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ABSTRACT 

Diffusive shock acceleration at collisionless shocks is thought to be the source of many of the ener- 
getic particles observed in space. Large-scale spatial variations of the magnetic field has been shown 
to be important in understanding observations. The effects are complex, so here we consider a sim- 
ple, illustrative model. Here, we solve numerically the Parker transport equation for a shock in the 
presence of large-scale sinusoidal magnetic-field variations. We demonstrate that the familiar planar- 
shock results can be significantly altered as a consequence of large-scale, meandering magnetic lines of 
force. Because perpendicular diffusion coefficient n± is generally much smaller than parallel diffusion 
coefficient km , the energetic charged particles arc trapped and preferentially accelerated along the 
shock front in the regions where the connection points of magnetic field lines intersecting the shock 
surface converge, and thus create the "hot spots" of the accelerated particles. For the regions where 
the connection points separate from each other, the acceleration to high energies will be suppressed. 
Further, the particles diffuse away from the "hot spot" regions and modify the spectra of downstream 
particle distribution. These features are qualitatively similar to the recent Voyager's observation in the 
Heliosheath. These results are potentially important for particle acceleration at shocks propagating in 
turbulent magnetized plasmas as well as those which contain large-scale nonplanar structures. Exam- 
ples include anomalous cosmic rays accelerated by the solar wind termination shock, energetic particles 
observed in propagating hcliospheric shocks, and galactic cosmic rays accelerated by supernova blast 
waves, etc. 

Subject headings: acceleration of particles - cosmic rays - shock waves - magnetic field 



1. INTRODUCTION 

Collisionless shocks in space and other astrophysi- 
cal environments are efficient accelerators of energetic 
charged-parti cles. Diffusive s hock acceleration ( here- 
inafte r PSA: IKrvmskvl [l977t lAxford et all fl977t iBelil 
119781: iBlandford fc Ostrikeri 119781) . is the most popular 
theory for charged-particle acceleration. It naturally pre- 
dicts a universal power-law distribution / oc p -7 with 
7 ~ 4.0 for strong shocks, where / is the phase-space 
distribution function, close to what observed in cosmic 
rays in many different regions of space. The basic con- 
clusions o f PSA can be drawn from the Parker transport 
equation (lParkerlll965D by considering the shock to be a 
compressive discontinuity in an infinite one-dimensional 
and time steady system. PSA is thought to be the mech- 
anism that accelerates anomalous cosmic rays (ACRs) in 
the Hcliospheric termination shock and also galactic cos- 
mic rays (GCRs) with energy up to at least 10 15 eV in 
supernova blast waves. However, recent in situ obser- 
vations in the termination shock and the Heliosheath 
by Voyager 1 (|Stone et al.l 120051 ) found the intensity 
of ACRs is not peaked at the termination shock and 
the energy spectrum is still unfolding after entering the 
Heliosheath, which strongly indicates the simple planar 
shock model is inadequate to interpret the acceleration 
of ACRs. Numerical and analytical studies suggest the 
possible solution can be made by considering the tern 



porary and/or spatial va r iation dFlorinski fe Zankl 



McComas fc Schwadronl I2006L Tjokipii fc Kotal 



2006 



2008 



Kota fc Jokipiil 120081: ISchwadron et al.ll2008D . In partic 
ular. iMcComas fc Schwadronl (|2006l ) discussed the im- 



portance of the magnetic geometry of a blunt shock 
on particle acceleration. They argued that the miss- 
ing ACRs at the nose of the Hcliospheric termination 
shock is due to particle energization occuring primar- 
ily back along the flanks of the shock where magnetic 
field lines have had a longer connect i on tim e and higher 
injection efficiency. iKota fc Jok ipii (2008|) presented a 
more sophisticated sim ulation which gives resul t s sim i- 
lar to that descri b ed by IMcComas fc Schwadronl (|2006|) . 
ISchwadron et al.l (|2008l ) also developed a 3-P analytic 
model for particle acceleration in a blunt shock, includ- 
ing perpendicular diffusion and drift motion due to large- 
scale shock structure. 

Large-scale magnetic field line meandering is ubiqui- 
tous in th e heliosphere and other astrophysical envi- 
ronmcnts (j Jokipiil [19661: Uokipii fc ParW[l96l iParkerl 
1979). The acceleration of charged-particles in col- 
lisionless shocks has been shown to be strongly af- 
fected by m agnetic -fie ld turbulence at different scales 
dGiacalond l2005allbl: f Giacalone fc Neugebauerl [200l 
IGuo fc Giacalondl2010D . The large-scale magnetic field 
variation will have important effects on the shock accel- 
eration since the transport of charged particles is differ- 
ent in the direction parallel and perpen dicular to the 
magn etic field, as shown in early work (| Jokipiil I1982L 
Il987t ). The blunt sho cks and shocks with fluctuating 
front (|Li fc Zankl 12000) which have the similar geome- 
try, are also relevant to this problem. In this study we 
consider the effect of the large-scale spatial variation of 
magnetic field on PSA. 
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2. BASIC CONSIDERATIONS AND NUMERICAL 
MODEL 

The diffusive shock acceleration (DSA) can be studie d 
by solving the Parker transport equation (lParkerlll965l ). 
which describes the evolution of the quasi-isotropic dis- 
tribution function f(xi,p,t) of energetic particles with 
momentum p dependent on the position Xi and time t 
including effects of diffusion, convection, drift, accelera- 
tion and source particles: 



dt 
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Here the diffusion coefficient tensor, 



df 
d lnp 
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Ui is the 
The drift 
, where w 



convection velocity and Q is a local source, 
velocity is given by = (pcw/3q)V x (B/B 2 
is the velocity of the particle, c is the speed of light, and 
q is the electric charge o f the parti cle. 

iKota fc Jokipiil (|2008f) and iKotal (|2010l manuscript to 
be submitted) considered analytically a model in which 
the upstream magnetic field was a plane (say, x, y), with 
average direction in the y direction. The x-component 
of the magnetic field was composed of uniform sections 
(straight field lines) alternating in sign, which were pe- 
riodic in y. They find "hot spots" and spectral effects 
which illustrate the effect of an upstream meandering in 
the magnetic field. 

Here, we consider a 2-D (x, z) system with a pla- 
nar shock at x = 0, and a sinusoid magnetic field 
B = Bo + sm(kz)5~B. For most of the parts in this 
paper we discuss the case shown in Figure 1. In this 
figure, the magnetic lines of force are illustrated by blue 
lines. The shock is denoted by the red dashed line. The 
system is periodic in the z direction, with the magnetic 
field convecting from upstream (x < 0) to downstream 
(x > 0). In the shock frame, the particles will be sub- 
jected to convection and diffusion due to the flow ve- 
locities U\ (upstream) and [^(downstream) and diffusion 
coefficients parallel and perpendicular to large scale mag- 
netic field («i(|i j_) and «2(||,_L))) respectively. The gra- 
dient and curvature drifts in this case are only in the 
direction out of the x — z plane and thus irrelevant to 
this study. Because of the steady velocity difference 
between upstream and downstream, charged particles 
which travel through the shock layer will be accelerated. 
However, since we consider the large-scale magnetic field 
variation, transport of energetic particles in the fluctu- 
ating magnetic field become important. The diffusion 
coefficient in the x — z system can be expressed as: 



Jij - 



(k_l - K\\)BjBj 
B 2 



(2) 



The normalization units chosen in this study are: the 
spatial scale Xq = 10 AU, the upstream velocity U\ = 
500 km/s, the time scale To = 3 x 10 6 sec and the dif- 
fusion coefficients are in unit of reo = 7.5 x 10 21 cm 2 /s. 
The shock compression ratio r = U1/U2 is taken to be 
4.0. The shock layer is considered to be a sharp vari- 
ation U x = (Ui + U 2 )/2 - (Ui - U 2 )tanh(x/th)/2 with 
thickness th = 1 x 10~ 3 , which is required to be less than 
Kxxx/Ui everywhere in the upstream simulation domain, 
where K X xl is the upstream diffusion coefficient normal 



to the shock surface. The simulation domain is taken 
to be [-2.0 < x < 2.0, —it < z < it]. The parallel 
diffusion coefficients upstream and downstream are as- 
sumed to be the same and taken to be k\\i = rep = 0.1 
at p = pq. The ratio between parallel diffusion coeffi- 
cient and perpendicular diffusion coefficient is taken to 
be k_l/km = 0.05, which is consistent with that deter- 
mined by integrating the tra jectories of test particles i n 
magnetic turbulence models (jGiacalone fc Joki"piilll999|) . 
The momentum dependence of the diffusion coefficient is 
taken to be re oc p 4 / 3 , corresponding to non-rel ativistic 
parti cles in a Kolmogorov turbulence spectrum () Jokipiil 
119711) . The time step is 1 x 10 -7 , which is small enough to 
resolve the variation of U x at the finite shock layer. We 
use a stochastic integration method, which is described 
in the Appendix, to obtain the numerical solution of the 
transport equation. The pseudo-particles are injected at 
the shock with initial momentum pq, and will be accel- 
erated if they across the shock. Particles which move 
past the upstream or downstream boundaries will be re- 
moved from the simulation. The system is periodic in the 
z— direction, so a pseudo-particle crossing the boundaries 
in z will rc-appcar at the opposite boundary and continue 
to be followed. A particle splitting technique similar to 
(|Giaca lone 2005al) is used in order to improve the statis- 
tics. Although we use very approximate parameters, we 
note that the results are insensitive to the precise num- 
bers. 

3. RESULTS AND DISCUSSION 

3.1. A shock propagating perpendicular to the average 
magnetic field 

Consider first the case where the average magnetic field 
is in the z-direction and the fluctuating magnetic field is 
SB = Bq. As shown in Figure 1, the magnetic field 
is convected through the shock front and is compressed 
in the x direction, thus B z2 = rB z \. For the sinusoid 
magnetic field considered in this paper, the local angle 
between upstream magnetic field and shock normal, Osn, 
will vary along the shock surface. As a magnetic field line 
passes through the shock surface, its connection points 
(the points where the field lines intersect the surface of 
the shock) will be moving apart in the middle of the 
plane (z = 0) and approaching each other on the both 
sides of the system (\z\ = n). Since, K|| 3> re_i_, the parti- 
cles tend to remain on the magnetic field lines. Because 
the acceleration only occurs at the shock front, as the 
magnetic lines of force convect downstream, the parti- 
cles will be trapped and accelerated at places where the 
connection points converge toward each other, leading to 
further acceleration. For the regions where the field lines 
separate from each other, the particles are swept away 
from those regions. Figure 2 displays the spatial distri- 
bution contours of accelerated particles in three energy 
ranges: 3.0 < p/p a < 4.0 (top), 8.0 < p/p a < 10.0 (mid- 
dle), and 15.0 < p/po < 30.0 (bottom). The density is 
represented by the number of particles in simulation and 
its unit is arbitrary. It can be seen that "hot spots" form 
in the regions that connection points approaching each 
other at all energy ranges, with lobes extend along the 
magnetic field lines. The density of the accelerated par- 
ticles at the connection-point separating region (in the 
middle of the plane) is clearly much smaller, although 
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there is still a concentration of low-energy accelerated 
particles there since the acceleration of low-energy par- 
ticles is rapid and efficient at perpendicular shocks. At 
higher energy ranges (middle and bottom), the lack of 
accelerated particles may be interpreted as due to the 
fact that the acceleration to high energies takes time. 

Figure 3 illustrates the profiles of the density of accel- 
erated particles for different energy ranges at z = (top) 
and z = 7r (bottom). In each panel, the black solid lines 
show the density of low energy particles (3.0 < p/po < 
4.0), the blue dashed lines show the density of interme- 
diate energy particles (8.0 < p/po < 10.0), and red dot 
dashed lines show the density of particles with high en- 
ergies (15.0 < p/po < 30.0). In connection-point sep- 
arating regions z — (top), it can be seen that while 
the downstream distribution of low energy particles is 
roughly a constant, the density of particles with higher 
energies increase as a function of distance downstream 
from the shock. These particles are not accelerated at 
the shock layer in the center of the plane but in the "hot 
spots". At z = 7r (bottom), the density of particles of 
all energies decreases as a function of distance, which in- 
dicates the accelerated particles diffusive away from the 
"hot spots" . Since the high energy particles have larger 
diffusion coefficients than the particles with low energy, it 
is easier for them to transport to the middle of the plane. 
The profile at z = is similar to Voyager's observation of 
anomalous cosmic ra ys (ACRs) at the termination shock 
and t he Heliosheath (jStone c t aT 200jl ICummings et al.l 
120081 ) which shows the intensity of the ACRs is still in- 
creasing and the energy spectrum is unfolding over a 
large distance after enteri ng the Helioshe a th. T he same 
physics has discussed by Uokipii fc Kotal (|2008t ). where 
"hot spot" of energetic particles is produced by the spa- 
tial variation of the injection of the source particles. In 
our current work, the concentration of energetic parti- 
cles are a consequence of particle accelerated in a shock 
containing large-scale magnetic variation. 

The top panel in Figure 4 represents the positions in z 
direction and the times as soon as the particles reached 
a certain momentum p c = 3po. We show that particles 
are accelerated mainly at the connection-point converg- 
ing region. There are also particles accelerated at middle 
of the plane because the particles can gain energy rapidly 
at perpendicular or highly oblique shock d ue to the smal l- 
ness of perpendicular diffusion coefficient ()Jokipiilll987t ). 
however, the further acceleration is suppressed by the 
effect that the charged particles travel away from the 
connection-point separating region, see also the top panel 
in Figure 5. It is clear that since the particles tend to 
follow the magnetic field lines, when the field line con- 
nection points separate from each other as field convects 
through the shock, the particles travel mainly along mag- 
netic field and away from the middle of the plane. The 
characteristic time for a field line convect from upstream 
to downstream r c = D/U\ ~ 1, therefore there is no 
significant acceleration in the middle of the plane after 
t = t c . Some of the particles can get more accelera- 
tion traveling from other region to "hot spots". Figure 
4 bottom shows the distance particles traveled in z di- 
rection from its original places \z — z$\ versus time when 
particle get accelerated at a certain energy (p c — 3po). It 
shows many particles are accelerated close to their orig- 
inal position, which is related to the acceleration in the 



"hot spot" . Nevertheless, there are also a number of par- 
ticles travel from the connection point separating region 
to "hot spot" and get further acceleration, which is rep- 
resented by the particles that travel a large distance in z 
direction. 

Figure 5 shows the same plot as Figure 4, except here 
the critical momentum is p c = lO.Opo- It is shown again 
in Figure 5 top that most of particles are accelerated to 
high energy are in the hot spot. However, as opposite to 
Figure 4 top, there are very few particles accelerated at 
the center of plane since energetic particles more trans- 
port away from the middle region and the time available 
is not long enough. For a quick estimate, the acceleration 
time is approximately, 

race = ^ ^(^f + ^f)rflnp 

Therefore for most of particles, they do not have suffi- 
cient time to be accelerated to high energies at the center. 
A number of particles accelerated at the center will travel 
to the "hot spot" and get more acceleration, as shown in 
Figure 5 bottom. 

Clearly, this presents different pictures of particle ac- 
celeration by the shock contain 2-D spatial magnetic 
field variations, indicates the resulting distribution func- 
tion is spatially dependent. In Figure 6 we show the 
steady state energy spectra obtained in the regions top: 
[0.1 < x < 0.3, 7r - 0.2 < z < it] (black solid line), and 
[0.1 < x < 0.3, -0.1 < z < 0.1] (green dashed line) and 
bottom: [0.8 < x < 1.0, n - 0.2 < z < n] (black solid 
line), and [0.8 < x < 1.0, -0.1 < z < 0.1] (green dashed 
line). It is shown that the spatial difference among dis- 
tribution functions at different locations caused by large- 
scale magnetic field variation is considerable. The black 
lines in both top plot and bottom plot, which correspond 
to the "hot spots" , show power-law like distributions ex- 
cept at high energies. At high energies, the particles 
will leave the simulation domain before gain enough en- 
ergy which causes the roll over in distribution function, 
this roll over is mainly caused by a finite distance to up- 
stream boundary. For other locations, the 2-D effect we 
discussed will produce the modification in distribution 
functions. The most pronounced effect can be found at 
the nose of the shock (green lines) , in the top panel the 
distribution of which shows a suppression of acceleration 
at all the energies. This insufficient acceleration is most 
prominent in the range of 6 — 12po- At these energies 
the acceleration time scales are longer than the time for 
the field line convection swipe the particles away from 
the connection-point separating region, as we discussed 
above. The bottom plot show that deep downstream the 
spectrum of accelerated particles is similar at high ener- 
gies since the mobility of these particles. 

3.2. An oblique shock 

The previous discussion has established the effect of a 
spatially varying upstream magnetic field on the accel- 
eration of fast charged particles at a shock which prop- 
agating normal to the average upstream magnetic field. 
We next consider the case where the shock propagation 
direction is not normal to the magnetic field. 
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Clearly, if the varying direction of the upstream mag- 
netic field is such that at some places the local angle 
of the magnetic field relative to the average field direc- 
tion exceeds the angle of the average magnetic field to 
the shock plane, we will have situations similar to that 
discussed in the previous sections. There will be places 
where the connection points of the magnetic field to the 
shock move further apart or closer together. Hence we 
expect the same physics can be applicable. An example is 
given in Figure 7. In this case the ratio of SB /Bo is taken 
to be 0.5, the averaged shock normal angle Osn = 70°. 
It can be seen from this plot that the connection points 
can still move toward each other in some regions. Figure 
8 shows the density contours of accelerated particles the 
same as Figure 2, but for the case of the oblique shock. 
We find that in this case the process we discussed in the 
last section is still persistent, even for a oblique shock 
and relative smaller SB/Bq. The "hot spot" forms corre- 
spond to the converging magnetic connection points and 
particle acceleration is suppressed in the region where 
connection points separate from each other. We may 
conclude that, for a shock which is oblique, if some mag- 
netic field lines can intersect the shock multiple times, 
we have "hot spots" of accelerated particles forms where 
the connection points converging together. 

4. SUMMARY AND CONCLUSIONS 

The acceleration of charged particles in shock waves 
is one of the most important unsolved problems in space 
physics and astrophysics. The charged particle transport 
in turbulent magnetic field and acceleration in shock re- 
gion are two inseparable problems. In this paper we illus- 
trate the effect of a large-scale sinusoidal magnetic field 
variation. This simple model allows a detailed examina- 
tion of the physical effects. As the magnetic field lines 
pass through the shock, the connection points between 
field lines on the shock surface will move accordingly. We 
find that the region where connection points approaching 
each other will trap and preferentially accelerate parti- 
cles to high energies and form "hot spots" along the shock 
surface, somewhat in analogy to the "hot spots" postu- 
lated bv iJokipii fc Kotal (|2008h . The shock acceleration 
will be suppressed at places where the connection points 
move apart each other. Some of the particles injected in 
those regions will transport to the "hot spots" and get 
further accelerated. The resulting distribution function 
is highly spatial dependent at the energies we studied, 
which could give a possible explanation to the Voyager 
observation of anomalous cosmic rays. Although we have 
discussed a simplified, illustrative model, the resulting 
spectra and radial distributions show qualitative simi- 
larity with the in situ Voyager observations. Thus, the 



intensities do not in general, peak at the the shock and 
the energy spectra are not power laws. We show this 
process is robust even in the case of oblique shocks with 
relatively small magnetic field variations. Large scale 
magnetic field variation, which could be due to magnetic 
structures like magnetic clouds, or the ubiquitous large 
scale field line random walk, will strongly modify the 
simple planar shock solution. This effect could work in a 
number of situations for large scale shock acceleration in- 
cluding magnetic variations, for example, the solar wind 
termination shock and supernova blast waves. 
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APPENDIX: STUDY DIFFUSIVE SHOCK 
ACCELERATION USING STOCHASTIC 
INTEGRATION METHOD 

In a 2-D system we considered, the Parker transport 
equation (1) can be re-written as, 



df 
dt 



df 

XX Q x 



Vxz dz 



dz \ xz dx 



as 



-U si 



1 dU x 
3 dx 



df 



dlnp 



(3) 



The equivalent stochastic differential equation is, 
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Where n, r 2 , and r 3 are different sets of random num- 
bers which satisfy < rj >= and < rf >= 1. In order to 
study diffusive shock acceleration, we approximate the 
shock layer as a sharp variation U x = (Z7i + 11%) / 2 — 
(Ui — U2)tanh(x/th)/2 with a thickness th much smaller 
compared with the characteristic length of diffusion ac- 
celeration Kxxi/Ui- At the same time, we have to make 
sure the time step At is small enough to resolve the mo- 
tion in the shock layer. 
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Fig. 1. — The shock and the magnetic field geometry for an upstream average magnetic field perpendicular to the shock normal. The 
blue lines represent the magnetic field lines and red dashed line indicates the surface of shock wave. The flow velocities are Ui (upstream) 
and U% (downstream). 
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Fig. 2. — The representation of density contour of accelerated particles, for energy range: 3.0 < p/po < 4.0(top), 8.0 < p/po < 
I0.0(middle), 15.0 < p/po < 30.0 (bottom). It is shown the hot spots forming on the both side of the system. The acceleration at the center 
of the shock is suppressed. 
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Fig. 3. — The profiles of density of the accelerated particles, for energy ranges: 3.0 < p/po < 4.0(black lines), 8.0 < p/po < 10.0(bluc 
dashed lines), 15.0 < p/po < 30.0 (red dot dashed lines) at different locations z = 0.0 (top) and it (bottom), respectively. 
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Fig. 5. — top: The position in z direction and time when the particle momentum reached p ■ 
direction \z — zq\ and time when the particle momentum reached p = lO.Opo- 



lO.Opo; bottom: The travel distance in z 
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Fig. 6. — The steady state distribution functions at Top: [0.1 < x < 0.3, n-0.2 < z < n] (black solid line), [0.1 < x < 0.3, -0.1 < z < 0.1] 
(green dashed line) and Bottom: [0.8 < x < 1.0, n — 0.2 < z < n] (black solid line), [0.8 < x < 1.0, —0.1 < z < 0.1] (green dashed line), 
respectively. 
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Fig. 7. — The shock and the magnetic field geometry considered for the case of an average magnetic field is 70° of the shock normal. The 
blue lines represent the magnetic field lines and red dashed line indicates the position of shock wave. The flow velocities are U\ (upstream) 
and t/ 2 (downstream). 



Shock Acceleration 



13 



3.0<p/p <4.0 




-2-10 1 2 
8.0<p/p <10.0 




-2-10 1 

15.0<p/p <30.0 



1.5506 
0.3636 
0.0853 
0.0200 

0.0047 
0.001 1 



■ 



0.0548 
0.0205 
0.0077 
0.0029 

0.001 1 
0.0004 



u 




Fig. 8. — The density contour of accelerated particles in the case of oblique shock and SB/Bo = 0.5, for energy range: 3.0 < p/po < 
4.0(top), 8.0 < p/p < 10.0(middle), 15.0 < p/p < 30.0{bottom). 



